Interpretation of the Flavor Dependence of Nucleon Form Factors 
in a Generalized Parton Distribution Model 



J. Osvaldo Gonzalez-Hernandez, Simonetta Liuti, 1 ^ 
Gary R. Goldstein, 2 '^ and Kunal Kathuria 1 '^ 

1 Department of Physics, University of Virginia, Charlottesville, VA 22904, USA. 
2 Department of Physics and Astronomy, 
Tufts University, Medford, MA 02155 USA. 

Abstract 

We give an interpretation of the u and d quarks contributions to the nucleon electromagnetic 
form factors for values of the four-momentum transfer up to few GeV where flavor separated data 
have been recently made available. The data show, in particular, a suppression of d quarks with 
respect to u quarks at large momentum transfer. This trend can be explained using generalized 
parton distributions which provide a correlation between momentum and coordinate spaces, both 
of which are necessary in order to interpret the partonic substructure of the form factors. We show 
that a flavor dependence originates from both Regge and quark-diquark terms in our model. 
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A recent experimental analysis displays a flavor separation of the nucleon elastic electro- 
magnetic form factors [lj. The data show both a suppression of the d quark contribution 
to the form factor with respect to the u quark at large momentum transfer. In this Letter 
we propose an interpretation of this somewhat puzzling behavior while simultaneously ex- 
ploring the correlation between momentum and coordinate space of the Generalized Parton 
Distributions (GPDs) or, more precisely, the connection they provide between single-particle 
densities in both momentum and impact parameter space. This approach enables us to show 
how the flavor dependence originates from different components of our model, namely a con- 
tribution from reggeon exchanges, and a handbag contribution, or quark-diquark term. 

GPDs define hybrid properties of the deeply virtual structure of nucleons [2H1] . Depend- 
ing on the ranges of the defining kinematical variables £ (skewness) and x (longitudinal 
momentum fraction), the GPDs, H(x,£,t), can be interpreted at a given scale, Q 2 , as 
non-forward Parton Distribution Functions (PDFs), x >| £ |, or Distribution Amplitudes 
(DAs), x <| £ |. Lorentz invariance implies that their first moment in the initial parton's 
longitudinal momentum fraction, x, is £ independent. It is the nucleon form factor, 



By Fourier transforming H(x,0,t) in the variable Aj_ (t = — Aj_), one obtains a single- 
particle density, or a diagonal object in impact parameter space, p(x,b). 

The formal backbone to this picture - the dominance of the handbag diagram - is provided 
by well established factorization theorems [3j|5]. The models which render a QCD description 
on how to interpret the parton correlators both in the forward and non-forward cases, and 
on how to establish the dominant degrees of freedom and their relative roles, are, however, 
debatable. Reviews on both the history and the more recent important developments on 
this longstanding issue, focusing on elastic scattering, can be found in e.g. [SHE]. 

Here we proceed by introducing a GPD perspective in this debate. In ^\ we developed a 
model in order to interpret Deeply Virtual Comption Scattering (DVCS) data. The model 
is a realization of the reggeized quark-diquark picture [TU] , in which Regge behavior of the 
GPDs is incorporated by introducing a spectral distribution for the spectator diquark mass. 
Upon integration over the mass, the spectral distribution yields on one side the desired x~ a 
behavior, and on the other, it is consistent with the diquark model (FigjT^t,b). 

The need for introducing a Regge term while applying diquark models to GPDs was 
realized in previous studies [IT] . There it was noticed that while it is a known fact that 
the diquark model cannot produce a steep enough increase of the structure functions at low 
x, this might be of minor importance in kinematical regions centered at relatively large x 
where most data in the multi-GeV region are. It is however a necessary contribution in order 
to obtain the normalization of the structure functions correctly. This observation becomes 
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important for GPDs where we require them to reproduce the form factor's behavior exactly 
through their normalization, or first moment as given in Eq.Q. The Regge term is therefore 
an essential ingredient in model building. The importance of the Regge term was realized 
recently also in Refs.[I21[I2]. There, however, the more singular behavior is introduced with 
a slightly different motive - it is required in order to model GPDs from a single Double 
Distribution. 

In practical terms, our reggeized diquark model depends on a number of parameters 
that we divide into Regge and pure diquark contributions, as we will explain in more detail 
below. The parameters were fixed by a fit applied recursively to PDFs from deep inelastic 
scattering data, and to form factors and DVCS data from Jefferson Lab |14j . The model 
was subsequently compared to data on different observables (charge and transverse single 
spin asymetries), in a different kinematical regime from HERMES [15j HB]- We define our 
parametrization as "flexible" in that, mostly owing to its recursive feature, the different 
components can be efficiently fitted separately as new data come in. 

Summarizing, we consider Ref. [9] as the accomplishment of a first phase in which we 
constructed a reggeized diquark model which satisfies fundamental requirements, such as 
polynomiality, positivity, crossing symmetries, hermiticity, and time reversal invariance. In 
the process, we studied the behavior of the various parameters both for the forward limit and 
for the integral relations including form factors, and we reproduced a number of observables. 
Our main result is summarized in Table I of Ref. [H] where an optimal set of the parameters 
obtained from data available at the time of publication was presented. 

With a viable model in hand we can now move on to understanding the role and interplay 
of its different components in interpreting a variety of experimental data. In particular, since 
the t dependence arises naturally in our model, namely it is not superimposed ad hoc, we 
can understand what features of the GPDs can simultaneously fit the PDFs and feed into 
the form factors. Why can, for instance, our flexible model reproduce the flavor dependence 
of the nucleon form factors? What components dominate - Regge or diquark, and for what 
values of t? What aspect of the nucleon substructure can this be traced back to? 

In Ref. pQ the u and d quarks contributions to the nucleon form factors were extracted for 
the first time up to values of t in the few GeV region, using isospin symmetric decomposition 
formulae, 

F l{2) = 2F 1(2) + F l(2) ( 2a ) 

F i(2) = 2iT (2) + F( {2) (2b) 

where, as usual, F^, and F^ are the Dirac and Pauli u and d quarks contributions to the 
proton form factor. The data show the somewhat surprising results that at the largest t the 
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FIG. 1: (a) i-channel Reggeon exchange diagram; (b) u-channel diquark exchange. The box has 
mass Mx, with spectral distribution p(M x ) as described in the text; (c) the same as (a) and (b) 
exhibiting the parton content; (d) Regge cut with its parton content corresponding to a diquark 
correlation. 



contribution of the d quark is suppressed relative to u, for both Fi and F 2 . 1 

A re- fit of our model using the new data displayed in Figj2] produced, first of all, the im- 
portant result that we reduced the errors of our parametrized form for the GPD considerably. 
Our model can be summarized in this expression, 



F(X, c, t) = MG^JX, C, t) R a / (X, c, t) 



(3) 



where F 



H, E\ the functions Gu A . m and R% ,a ' parametrize respectively, the quark-diquark 



Mx,m """^ 

and Regge contributions; X and ( are related to x and ^ by, X = (a; + £)/(l + ^), £ = 
2(/(2 — (), however, here we are interested in the ( = £ = O^X = x limit. The diquark 
components read, 



M\,m q 



Af 



d 2 k< 



[m g + Mx) (m g + Mx) + k ± ■ kj 



1 - 

d 2 k ± 
1 — x 



x [M 2 q (x) - kl/(l ~ x)f{M 2 q (x) - k ± /(l - x)} 2 

-2M/A 2 ± \{m q + Mx) k± ■ A± - (m Q + Mx) k ± ■ Aj 



[M 2 (x) - k 2 J(l - x)] 2 [M 2 (x) - k ± /(l - x)Y 



(4a) 



(4b) 



1 It should be also noticed that both trends appear by extrapolating the widely used form factors 
parametrization by Kelly [T7] to large t. 
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FIG. 2: Left Panel: Proton form factors t F q (above) and k~ t F% (below) plotted vs. —t, as 
obtained from our parametrization using Eq.Q, by fitting the data preceding Ref.[lJ. Right Panel: 
same as Left, including the data from Ref.pQ. The parametrization's error bands were obtained 
without taking into account parameters correlations. Experimental data on both panels are from 
Ref.pQ. 

where q = u, d, M is in GeV 4 , kj_ = kj_ - (1 - x)A ± , M 2 q (x) = xM 2 - x/{\ - x)M q2 - M| 2 , 
M being the proton mass. 
The Regge term is given by 

R = A f x -l«+«'(*)t} j (5) 

where a!(x) = a'(l — x) p . 

Notice that while Eqs.Q and ^ are displayed in a compact form that is suitable for 
a parametrization, they also represent the final outcome of a model integrating significant 
physical information from both quark-diquark correlations, and Regge behavior. 

Eq.(|3]) corresponds to a generalization of the diquark picture in which the mass of the 
diquark, M q x is not fixed, but has a spectral distribution, pn(Mx) oc (M 1 



q \2 [a+a'(x)t-l] 



X i 



poo 

F(X,C,t) = N I dM%M^'*- 1] Gfr x>m {z,Q,t) 
Jm x 

POO 

= Afx- [a+a ' {x)t] / dzz^ a+a '^G^ A m (x,0,t) 



=z/M% 

Afx-^ a '^G%iJx,0,t). (6) 



where the superscript q is omitted for simplicity, and the variable M\ in m has been 
replaced with z = xM\ (more details will be given in future work). 
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FIG. 3: (color online) Combined x and t dependences of the Regge contribution in our model 
(Eq.([5]) and discussion in text). Upper panel: trajectories with x-dependent parameters from our 
model. The dot dashed line is a comparison with the softening of the slope first suggested in 
Ref . [201 121] . Lowe panel: Regge term, Eq.Q obtained using the trajectories shown in the upper 
panel. 

The Regge term has different interpretations in different kinematical regions: 
£ = 0, t = As first shown in Ref. pU] in a simple covariant model for PDFs, choosing a 
spectral distribution of the form pr(Mx) oc gives rise to the behavior x _ct (°) upon 

integration over all M\. In Ref. [9] we illustrated the "reggeization" procedure by considering 
the forward, spin independent GPD, H q (x,0,0) = ff( function of a continuum of 

diquark masses. Aside from overall constant factors, we obtained the expression for £ = 
displayed in Eqs.(|5]) and (|6j). 

C = 0, small t The result from Ref. (9] can be trivially extended to t < 0, noticing that the 
integrals are done at fixed t. So long as a + a'(x)t remains positive, i.e. for —t < a/a'(x), 
the physical interpretation remains unaltered because the behavior of the M\ spectrum does 
not change considerably. This is the region spanned in DVCS-type experiments, where the 
reaction's four-momentum transfer Q 2 w (few GeV 2 ), and — t < 1 GeV 2 . 

( = 0, larger t Larger, but not asymptotic, negative t can still be consistent with Regge 
behavior. The x-dependence of the slope parameter, a + a'(x)t, is such that it softens the 
fall off of the trajectory with negative t. Although this term was initially introduced driven 
by the need to guarantee more flexibility in parametrizations upon Fourier transformation in 
A± [HJ HU HE] , we are now able to interpret it physically: a!(x) effectively describes Regge 
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FIG. 4: (color online) Diquark vertex function contributions to the Dirac and Pauli form factors 
described in text. 

behavior, including an important deviation from naive Regge, namely the appearance of cuts. 
The behavior with t is illustrated in Fig|3j It is interesting to relate this to the partonic 
structure that gives rise to Regge behavior represented in FigJT]3,d. Fig. corresponds to 
a simple Reggeon exchange while Fig. [l]i exhibits a Regge cut. Notice that the non-planar 
structure of this graph allows us to interpret the intermediate region as given by coupling 
of diquarks to the virtual photon. A similar contribution was considered in the Poincare 
covariant Dyson-Schwinger equation (DSE) approach of Ref. [6], which includes both dressed 
quark and diquark components coupling to the virtual photon. 

( 7^ Including skewness, ( ^ is more complicated, as singularities at the crossover 
point, X = ((x = £), might arise. This issue has been extensively discussed recently in [12] 
in the context of Double Distributions (see also [13]), although in the t = limit. Detailed 
evaluations of the ( ^ case do not enter the form factors picture, and will be presented 
elsewhere. 

We now turn to the quark-diquark component by highlighting the angular momentum 
structure of the diquark correlations, given in this case by the J p = + (scalar) and J p = 1 + 
(axial vector) configurations. From these one obtains distinct predictions for the u and d 
quarks contributions using a SU(4) symmetric wave function for the proton (see e.g. [22]), 

F\X, C, t) = V=°(X, C, t) + V =1 (X, C, t) (7) 
F d (X,C,t) = F s=1 (X,C,t). (8) 
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From the practical point of view, one can take a similar form for the scalar and axial- 
vector couplings (scalar-like), while distinguishing their different contributions by varying 
their respective mass parameters as in Eqs.Q. This choice is motivated by previous work 
[TT] showing that the full account of the axial-vector coupling does not sensibly improve 
the shape of parametrizations, while considerably increasing the algebraic complexity of the 
various structures. 

Another important choice is in the form of the coupling at the proton-quark-diquark 
vertex. Consistently with studies of baryons in the context of the application of Dyson - 
Schwinger equations in QCD (|23j and references therein), we chose a dipole type coupling, 

k 2 -m 2 n 

r = g s |— , 9 

y (P-Mf ) 2 V ; 



where both m q and Mj[ are among our fitted parameters. In Fig|4] we show the integrated 
wave function contributions, J dx<p*(x, k — A)<j)(x, k), that enter the form factors calculation 
for a given value of t = — 2 GeV 2 . One can see that for the Dirac form factor the u quark 
tends to dominate over the d quark while for the Pauli form factor the u and d quark 
contributions are similar to one another. The t dependence of this term is crucial in our 
model: at low t one finds a harder d quark component than for the u quark (analyses along 
these lines were initiated in Refs. [2"H 25J, where, however, the flavor dependence was not 
addressed) . 

The behavior in the different regions of t discussed above, as well as the di quark contri- 
bution to the nucleon form factors can be observed in the data through our fitting procedure 
p. In a first step we fitted H q (x, 0, 0), q = u, d to the PDFs, q(x), and obtained the parame- 
ters in Eqs.Q and (jsj, Ml (Eq.([9|), M 9 X (the starting point of the diquark mass spectrum), 
m q (the quark mass), and a q . The fit was performed by evolving to the scale, Q 2 of the 
deep inelastic data (in most formulae herewith we omit displaying Q 2 explicitly because it 
is not relevant to the discussion of the form factor). We used the same set of parameters 
for E q (x, 0, 0) since this is unconstrained at t = (notice however that its functional form 
is different, Eq.(4b)). 

Keeping these parameters fixed we obtained a' and p by fitting the Dirac and Pauli form 
factors 

if B) (t) = J dxH^\x,0,t) = e u{d) F^ d \t) - e d{u) Ff {u \t), (10) 
if = / dxBF^ixM = e u(d) F^ d \t) - e d{u) Fi [u \t), (11) 

where e u = 2/3, and e d = —1/3. In Table [T] we show both sets of values for a' and p, using 
both the old set of data and the data from Ref. pQ. Results for both cases are displayed in 
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x 2 
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0.5 
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TABLE I: Parameters obtained from our recursive fitting procedure applied to H q , E q , q = u,d. 
We obtained a' q , p q , by fitting the proton and neutron electromagnetic form factors. Also shown 
are the \ 2 values for the separate contributions to the fit. 



1.2 
1 

0.8 

{ 0-6 
0.4 
0.2 

0.25 
0.2 
0.15 

V 

o.i 

0.05 




Regge 

I Total u 

T Total d 










U 






.. . 






f — 








u 






d 








1 i , , i 1 i , i , 1 , i i , 1 , , , , ] , i , i 1 , i , i 1 , , , i 1 i , i , 1 i , i , 



------ °" 



u. 

V 



0.5 1 1.5 2 2.5 3 3.5 4 4.5 
-t (GeV 2 ) 



1.2 
1 

0.8 
0.6 
0.4 
0.2 


1 

0.8 
0.6 
0.4 
0.2 




_ Diquark 


Diquark 


7 Regge 


Regge 


Total u 


Total d 


i \ 






^ \ 






- \ N. 


— 




iX 

"* \ 
\ 1 

— 1 \ 







1 



-t (GeV z ) 



-t (GeV") 



FIG. 5: (color online) Left Panel: i F', q = u,d (above) and k" rF 2 (below) plotted vs. — t, as 
obtained from our parametrization using Eq.([3]). Right Panel: Proton form factors (above), 
u and d quarks contributions displayed in the left and right panels, respectively; k" 1 ^ (below) 
Experimental data from Ref.jlj. The dotted lines correspond to the diquark contribution, Eq.(|13|); 



the dashed lines to the Regge contribution, Eq.(14). 



Figj2j 

The nucleon form factors are decomposed in terms of w-quark and <i-quark contributions, 
which in turn include a quark- diquark, and Regge term. The quark-diquark term describes 
scattering from a single, non-interacting u or d quark leaving either a ud diquark with spin 
S = 0, 1, or a or a uu diquark system with spin S = 1 as spectators. The Regge term 
describes scattering from a state which is effectively rendered coherent through Regge cuts, 
or diquark correlations (Fig|2]d) . 

In order to ascertain whether the Regge component alone can drive the t-dependence, or 
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FIG. 6: (color online) Left Panel: Integrated parton density distribution in transverse space for 
the Dirac form factor, displaying all components of our model. Right Panel: same as left, for the 
Pauli form factor. 

to determine its interplay with the t-dependent diquark term, we separated the Regge and 
diquark terms in the following way 

F q {X,C,t) = F£ + F! t (12) 

where 

n = (13) 

n = m +R RT' q (14) 

with f g(r) = h g{r) ,e g{r)] at g+r = (i/R a r'« + i/<|, m J _1 - 

Results for the Dirac and Pauli form factors calculated using Eqs. 

F? = fdx{H G + H R l (15) 
Jo 

Ft = [ dx(E G + E R ), (16) 

are shown in Fig|5j We notice that for the w-quark the Regge term does dominate the 
low t behavior as expected. For Ff the contributions of the Regge and diquark terms 
are instead comparable. The persistence of the Regge term at larger t where the form 
factors are dominated by large x [24J is not unexpected in our model, in view of diquark 
correlations/Regge cuts which extend the validity of the Regge model to this region (see 
Fig|3]and Refs. [201 EH)- 

The Fourier transforms of GPDs with respect to Ai define the parton density distribu- 
tions at a transverse position b for a given longitudinal momentum fraction, x, namely two 



Hp), 
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dimensional distributions in the transverse plane with respect to the proton's direction of 
motion. We can therefore connect each form factor component to hadronic distances from 
the proton's center of momentum defined as [27J, 
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(17) 
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d 2 A A 
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(b 2 )l G +{b 2 ) q 2R . 



2\q 



(18) 



In Fig|6| we show the density integrated over x, p1( 2 )(b) f° r & U components. The behavior 
of the form factors at larger t is reflected in the behavior of p as b — > 0. One can deduce 
that Dirac hadronic sizes for the u quarks (Fig|6j left) are smaller than for the d quarks. A 
smaller difference between u and d quarks appears for Pauli hadronic sizes. For the quark- 
diquark component of our model one can easily trace back this behavior to the transverse 
momentum distribution, Fig|3j which is harder for the u quark than for the d quark. This is 
in turn a consequence of the difference in masses between the axial vector and scalar diquark 
which respectively govern the d quark and u quark contributions. Similar conclusions were 
reached in Ref. [6] although from different symmetry properties since they include also the 
longitudinal spatial component. The Regge term in our model is less straightforward to 
interpret, however, if we view it as a coherent state, then the distance involved in this case 
would be the combination of the two-quark, composite, hadronic component's size containing 
either the struck u or d quark, and the overall distance of this component from the proton's 
center of momentum. 



In conclusion, we used a GPDs based approach as a way to understand the u and d quarks 
components of the nucleon form factors. The GPDs were evaluated using a reggeized quark- 
diquark model whose parameters are fixed to simultaneously fit the deep inelastic limit, the 
nucleon form factors, and DVCS data. In our approach both a Regge-type contribution 
and quark-diquark components determine the behavior of the form factors in the few GeV 
region. The Regge term is interpreted as effectively taking into account coherent processes, 
including final state interactions or Regge cuts. The latter are connected through duality 
to diquark correlations. As t increases we expect the Regge cuts to become suppressed, and 
the nucleon form factor to be dominated by Regge asymptotic behavior modulated by a t 
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dependent term. The quark-diquark term, on the contrary, involves scattering from a single 
quark, leaving a diquark system as a spectator (what happens in the center of the nucleon 
is determined by asymptotia which are not at reach in the kinematical domain of the data). 

A first outcome is that the new highly precise form factor data produce much improved 
constraints on our GPDs parameters. The interpretation of the flavor dependence of the 
data lies in the non-perturbative structure of both the Regge and quark-diquark terms. It is 
a subtle combination of effects that cannot be ascribed to a single, simply motivated mech- 
anism. We see on one side the persistence of notably flavor dependent Regge contributions 
at large t. Isovector exchanges present in our model are e.g. known to couple differently to 
u and d quarks. These terms alone produce a suppression of the d quark contribution with 
respect to the u quark. The Regge mechanism dominates the Pauli form factors. 

Similarly, scattering from a single quark through the quark-diquark mechanism is fla- 
vor dependent, the d quark being suppressed because of the softer transverse momentum 
dependence of the proton-quark-diquark vertex function. 

Through the concept of GPDs the Regge and diquark mechanisms are realized corre- 
spondingly in coordinate space, in the transverse plane where our description reflects the 
behavior of the form factors. Finally, a unified picture of the Regge and diquark contribu- 
tions can be given using duality arguments. The Regge terms include a component that 
corresponds to diquark correlations in the nucleon. In the GPD picture one can also study 
the internal spatial distribution and size of these components. 

The newly available flavor separated form factor data at large t stimulated this work. 
Further studies exploring a possibly important role of diquark/few-parton correlations inside 
the proton will be carried out in the near future. 

We thank Gordon Cates for helpful discussions. This work has been supported by the 
U.S. Department of Energy grants DE-FG02-01ER4120 (J.O.G.H., , K.K., S.L.), and DE- 
FG02-92ER40702 (G.R.G.). 
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